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Abstract — This brief paper further investigates the locally 
and globally adaptive synchronization of an uncertain complex 
dynamical network. Several network synchronization criteria 
are deduced. Especially, our hypotheses and designed adaptive 
controllers for network synchronization are rather simple in 
form. It is very useful for future practical engineering design. 
Moreover, numerical simulations are also given to show the 
effectiveness of our synchronization approaches. 

Index Terms — Complex networks, adaptive synchronization, 
uncertain systems 

I. Introduction 

OVER THE PAST DECADE, complex networks have 
been intensively studied in various disciplines, such as 
social, biological, mathematical, and engineering sciences [1- 
8]. A complex network is a large set of interconnected nodes, 
where the nodes and connections can be anything. Detailed 
examples are the World Wide Web, Internet, communication 
networks, metabolic systems, food webs, electrical power 
grids, and so on. 

Recently, one of the interesting and significant phenomena 
in complex dynamical networks is the synchronization of all 
dynamical nodes in a network. In fact, synchronization is 
a kind of typical collective behaviors and basic motions in 
nature. For example, the synchronization of coupled oscillators 
can explain well many natural phenomena. Furthermore, some 
synchronization phenomena are very useful in our daily life, 
such as the synchronous transfer of digital or analog signals 
in communication networks. Specifically, synchronization in 
networks of coupled chaotic systems has received a great deal 
of attention. Some synchronization criteria of two or three 
Lorenz systems have been obtained in the literature. However, 
it is often difficult to get the exact estimation of the coupling 
coefficients since we do not know the exact boundary for 
most chaotic systems. Up to now, we can only estimate the 
boundary of very few chaotic systems [9-13], such as the 
Lorenz, Chen, Lii systems [14]. Moreover, we often know 
very little information on the network structure, which makes 
network design very difficult. To overcome these difficulties, 
an effectively adaptive synchronization approach is proposed 
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based on an uncertain complex dynamical network model in 
this paper. 

Slotine, Wang, and Rifai [16,17] further discussed the 
synchronization of nonlinearly coupled continuous and hybrid 
oscillators networks by using the contraction analysis approach 
[18]. Bohacek and Jonckheere [19-20] proposed the so-called 
linear dynamically varying method based on discrete time 
dynamical systems. In the following, by using Lyapunov sta- 
bility theory, several novel locally and globally asymptotically 
stable network synchronization criteria are deduced for an 
uncertain complex dynamical network. Compared with some 
similar results [3,5,15], our sufficient conditions for network 
synchronization are rather broad and the controllers are very 
simple. It is very useful for future practical engineering design. 
Moreover, our analysis method and network model are very 
different from those of the above referenced literature [16-20]. 
However, for some complex systems (e.g., biological systems) 
with unknown couplings, our conditions are hard to be verified. 
In fact, it is impossible to propose a universal synchronization 
criterion for various complex networks since there are many 
uncertain factors, such as network structures and coupling 
mechanisms. 

This paper is organized as follows. An uncertain complex 
dynamical network model and several necessary hypotheses 
are given in Section II. In Section III, locally and glob- 
ally adaptive synchronization criteria for uncertain complex 
dynamical networks are proposed. In Section IV, a simple 
example is provided to verify the effectiveness of the proposed 
method. Finally, conclusions are given in Section V. 

II. Preliminaries 

This section introduces an uncertain complex dynamical 
network model and gives some preliminary definitions and 
hypotheses. 

A. An uncertain complex dynamical network model 

Consider an uncertain complex dynamical network con- 
sisting of N identical nonlinear oscillators with uncertain 
nonlinear diffusive couplings, which is described by 

ki = f(xj, t) + hi(xi, x 2 , • • • , x N ) + u i; (1) 

where 1 < i < N, Xi = (x il; x l2 , ••• , x ln ) T € R" is 
the state vector of the zth node, f : ft x R + — * R™ is a 
smooth nonlinear vector field, node dynamics is x = f(x, t), 
hi ; O X • • • X fi — > R" are unknown nonlinear smooth 
diffusive coupling functions, Ui G R" are the control inputs, 
and the coupling-control terms satisfy hj(s, s, • • • , s) + U, 
0, where s is a synchronous solution of the node system x = 
f(x, t). 
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B. Preliminaries 

Network synchronization is a typical collective behavior. In 
the following, a rigorous mathematical definition is introduced 
for the concept of network synchronization. 

Definition 1: Let Xj(i ; to, Xo) (1 < i < N) be a solution 
of the dynamical network {!}, where Xo = (x°, x 2) • • • , Xjy), 
f : n x R+ R", and h, : fi x • • • x fi -> R™ (1 < 
i < AT) are continuously differentiable, 57 C R™. If there is 
a nonempty subset A C f2, with x- e A (1 < j < N), such 
that Xi(t ; t , X ) G SI for alH > i > 1 < * < AT. an d 



lim [[Xi(* ; X ) - s(f ; to, x ) 



1 < % < N (2) 



where s(t; t n , x ) is a solution of the system x = f(x, t) 
with x G 51, then the dynamical network Q is said to 
realize synchronization and Ax ■ • • x A is called the region 
of synchrony for the dynamical network Q. 

Hereafter, denote s(t;to,x n ) as s(t). Then S(t) — 
(s T (t), s T (t), ■ ■ ■ , s T (t)) T is a synchronous solution of un- 
certain dynamical network Q since it is a diffusive coupling 
network. Here, s(t) can be an equilibrium point, a periodic 
orbit, an aperiodic orbit, or a chaotic orbit in the phase space. 

Define error vector 



e,(t) = Xi(t) - s(t), 1 < i < N. 



(3) 



Then the objective of controller Ui is to guide the dynamical 
network ([0 to synchronize. That is, 

lim \\e t {t)\\ 2 = 0, 1 < i < N. (4) 

t — > + oo 

Since s = f(s, t), from network Q, we have 

e !: = f(xj, s, t) + hj(xi, x 2 , • • • , xat, s) + Ui , (5) 
where 1 < i < A 7 ", f(x,;, s, t) = f(x,;, i) — 

f(s, t), &i(Xl, X 2 , • • • , X N , S) = hj(Xl, X 2 , • • • , Xat) - 

hj(s, s, ■ ■ • , s). 

In the following, we give several useful hypotheses. 

Hypothesis 1: (HI) Assume that there exists a nonnegative 
constant a satisfying ||Df(s, t)\\ 2 = ||A(t)|| 2 < a, where 
A(t) is the Jacobian of f(s, t). 

Hypothesis 2: (H2) Suppose that there exist non- 
negative constants 7^ (1 < i, j < N) satisfying 

N 

, xat, s) II 2 < E Tij ll e j II2 for 1 < i < N. 



|hi(xi, x 2 , • 
Remark 1: If HI holds, then we get 



A(t)+A T (t) 
2 



< a. 



III. Adaptive synchronization of an uncertain 

COMPLEX DYNAMICAL NETWORK 

This section discusses the local synchronization and global 
synchronization of the uncertain complex dynamical network 
0. Several network synchronization criteria are given. 



A. Local Synchronization 

Linearizing error system Q around zero gives 

ej = A(t)ei{t) + bi(xi, x 2 , • ■ • , x N , s) + Uj , (6) 

where 1 < i < N and recall that A(t) = Df(s, t) is the 
Jacobian of f evaluated at x = s(t). 

Based on HI and H2, a network synchronization criterion 
is deduced as follows. 

Theorem 1: Suppose that HI and H2 hold. Then the 
synchronous solution S(t) of uncertain dynamical network Q 
is locally asymptotically stable under the adaptive controllers 



-did, 1 < i < N 



(7) 



and updating laws 



kieja = hWag, 1 < i < n, (8) 

where fc, (1 < i < N) are positive constants. 
Proof: Define a Lyapunov candidate as follows: 



V 



1 N 

IE 



T 

e, e 4 



1 N 

-Y 

i=l 



(di - d t y 



(9) 



where di (1 < i < N) are positive constants to be deter- 
mined. Thus one gets 

N N - ■ 

V = I E (efe 4 + efe.) + £ 

i= 1 i= 1 

N 



Eef(^^-^I„)e i 

N 



N 



+ E ef h 4 (xi, x 2 , ••• , x N , s) + E (d; - ^i) e f e i 



< 



J- e T / A(t)+A T (t) 

! = 1 ' ^ 



- dil n e 



iV AT 

+ E E TiiNlallejHa 

i = 1 j = 1 
AT N N 

< E(«- di)||ei||l + E E 7^||ei||2||e i || 2 

i = l £ = 1 j = 1 

= e T (r + diag{a — di, a — d 2 , • • • , a — d N })e , 
where e = (||ei|| 2 , ||e 2 || 2 , • • ■ , ||eAr|| 2 ) T and T = 

(7ij)N x N- 

Since a and 7^ (1 < i, j < A 7 ) are nonnegative constants, 
one can select suitable constants di(l < i < N) to 
make T + diag{a — di, a — d 2 , • • • , a — 5n} a neg- 
ative definite matrix. Thus it follows that the error vector 
T] = (ej , ej, ■ • ■ , eJ5j) T — > as t — > +00. That is, the 
synchronous solution S(t) of uncertain dynamical network Q 
is locally asymptotically stable under the adaptive controllers 
and updating laws 0. 

The proof is thus completed. 

Assume that the coupling of network Q is linear satisfying 

jv 

hj(xi, x 2 , • • • , x^) = E bijXj for 1 < i < N, where 

3=1 

bij (1 < i, j < A 7 ) are constants. Then the uncertain network 
is recasted as follows: 



N 



±i = f(Xj, t) + ^ l '-J X ' + U i > 

3 = 1 



1 < i < N. 



(10) 
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For linear coupling, H2 is naturally satisfied. Thus one gets 
the following corollaries. 

Corollary 1: Suppose that HI holds. Then the synchronous 
solution S(t) of the uncertain dynamical network dlOt is 
locally asymptotically stable under the adaptive controllers Q 
and updating laws (|8jl. 

Moreover, for the coupling scheme hj(xi, X2, • • • , Xjv) 

Y 

E bijjf(Xj) with 1 < i < N, where b li (1 < i, j < N) 

N 

are constants satisfying E bij — for 1 < i < N and 

i=i 

j | Dp (^) 1 1 2 < S for £ S ft, the network l[T} is rewritten as 
follows: 

JV 

ii = t(xi, t) + fciPfo) + Ui) 1 < i < N. (11) 

i = i 



If HI holds, then one has 



JV 



||hi(xi, x 2 , ••• , Xjv, s)|| a - E IM IIp( x j) ~ P(s)|| 2 

i = i 

JV 

< E ^ lll^|| 2 

i = i 

for 1 < i < N. That is, H2 holds and one gets the following 
corollary. 

Corollary 2: Assume that HI holds. Then, under the adap- 
tive controllers and updating laws 0, the synchronous 
solution S(t) of the uncertain dynamical network dl 11 1 is 
locally asymptotically stable. 

In the following subsection, we discuss the global synchro- 
nization case. 



B. Global synchronization 

This section presents two global network synchronization 
criteria. 

Rewrite node dynamics x$ = f(xj, t) as Xj = Bxj(£) + 
g(xj, t), where B G R nX " is a constant matrix and g : 
n x R+ -» R" is a smooth nonlinear function. Thus network 
Q is described by 

Xi = Bxi(t) + g(xi, t) + hj(xi, x 2 , • • • , x N ) + Ui, (12) 

where 1 < i < N. Similarly, one can get the error system 

ei = Bei(t) + g(Xi, s, t) + hj(xi, x 2 , • • ■ , xat, s) + u 4 , 

(13) 

where 1 < i < TV and g(x i5 s, t) = g(x», i) - g(s, i). 

Hypothesis 3: (H3) Suppose that there exists a nonnegative 
constant fi satisfying ||g(x», s, t)\\ 2 < Mlkilb- 

Then one can get the following global network synchroniza- 
tion criterion. 

Theorem 2: Suppose that H2 and H3 hold. Then the 
synchronous solution S(t) of uncertain dynamical network Q 
is globally asymptotically stable under the adaptive controllers 



and updating laws 



&i 



Nil- 



1 < i < N, 



(15) 



where ki (1 < i < AT) are positive constants. 

Proof: Since B is a given constant matrix, there exists a 
nonnegative constant (3 such that ||B||2 < (3. It follows that 

< 0- 



B + B J 



Similarly, construct Lyapunov function (|9j, then one has 



V 



N 
i = 1 



+ B 1 



N 



E efgj(Xi, s, t) 



E ef hi(xi, x 2 , 

i = l 



XiV, s) 



< 



IV 

ECS 

1 = 1 
e T (T h 



h M - djJIIeiHi + E E TtfNWMa 

i=lj=l 

diag{/3 + - di, • • • , (3 + fi - d N })e . 



ei 2, e 2 2 , 



\e N \\ 2 ) T and T 



where e = 
(7u)n x N- 

Since f3, fj, and 7^ (1 < i, j < N) are nonnegative 
constants, one can select suitable constants di (1 < i < N) 
to make T + diag{/3 + fj, — di, f3 + fi — d 2 , • • • , /3 + 
/i — 5n} a negative definite matrix. Then the error vector 

T „T T\T 



i] 



(e 



1 I c 2 I 



as t 



-00. That is, the 



synchronous solution S(t) of uncertain dynamical network 
is globally asymptotically stable under the adaptive controllers 
H41 and updating laws Jl 5b . 
This completes the proof. 

Similarly, one gets the following two corollaries of global 
network synchronization. 

Corollary 3: Suppose that H3 holds. Then the synchronous 
solution S(i) of uncertain linearly coupled dynamical network 
i ll Oi l is globally asymptotically stable under the adaptive con- 
trollers dl4> and updating laws J15I . 

Corollary 4: Suppose that H3 holds. Then the synchronous 
solution S(t) of uncertain dynamical network dl Q is globally 
asymptotically stable under the adaptive controllers dl4t and 
updating laws dl5t . 

Proof. According to dl Q , one has 



|hi(xi,x 2 , • ■ • ,xjv,s)|| 2 



N 

E bijiBej 



< 



E b 



ij 1 



ilBll 



gj(Xi,S,t)) 



-did, 1 < i < N 



(14) 



thus H2 holds. Therefore, from Theorem 2, the synchronous 
solution S(t) of network dl Q is globally asymptotically stable 
under the adaptive controllers J14b and updating laws dl5t . 
The proof is thus completed. 

Hypothesis 4: (H4) Assume that g(x, i) satisfies the Lip- 
schitz condition. That is, there exists a positive constant k 
satisfying ||g(x, t) — g(y, t)\\ < k\\x — y\\, where k is the 
Lipschitz constant. 

Obviously, H4 implies H3. Now one has the following 
synchronization criterion. 

Theorem 3: Suppose that H2 and H4 hold. Then the 
synchronous solution S(t) of uncertain dynamical network Q 
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is globally asymptotically stable under the adaptive controllers 
dl4> and updating laws dl5l . 



IV. An example 

This section presents an example to show the effectiveness 
of the above synchronization criteria. 

Consider a dynamical network consisting of 50 identical 
Lorenz systems. Here, node dynamics is described by 



7 l] 








' ° 1 






Xi2 


+ 


%il%i3 1 


V %i3 j 




\ x l3 ) 







where 







a = 10, b = | , c = 28 and 1 < i < 50. And the 
networked system is defined as follows: 

xn \ I 

A I Xi2 + -XilXt3 
Xi 3 J \ XiiX i2 

/i(xi) - 2/x(x i+ i) + /i(x l+2 ) 


/ 2 (x. i )-2/ 2 (x l+1 ) + / 2 (x. i+2 ) 



and 



dj — fcj ||Cj||2 j 



(16) 
(17) 



= a(x i2 - x a ), / 2 (Xj) = x a x i2 - bx l3 , x 51 
Xi, x§2 = x 2j an d 1 < * < 50. 
Obviously, one gets 



g(Xi, s, t) = 





~Xi\X i3 + SiS 3 
Xi\Xi2 — S\S2 



-Xi 3 en - sie a 
Xi 2 en + sie i2 



where 1 < i < 50. 

Since Lorenz attractor is confined to a bounded 
region $ C R 3 [9-13], there exists a constant 
M satisfying \x. 



< M for 1 < 



50 



an d j = 1, 2, 3. Therefore, [| g(x f , s, t)|| 2 
y/Jxaeu + sie i3 ) 2 + (x i 2e l i + sie l2 ) 2 < 2Af||ei|| 2 . 
Similarly, one has 



hi(xi,x 2 , • • • ,xjv,s) 



/i(ei) - 2/i(e i+ i) + /i(e i+2 ) 


/3( x ii x i+l, x i+2, S) 



-be 



i3 



2be 



t+1,3 



be 



i+2,3 



where / 3 (x i; x i+1 , x J+2 , s) 
Xi2Bn + S\ei2 — 2(x i+ i , 2 e i+11 + Sie^+i^) + x i+2 ,2^i+2,i 
siei+2.2 and 1 < i < 50. 




1 2 

(a) e a (1 < i < 50) 




1 2 

(b) e l2 (1 < i < 50) 




(c) e i3 (1 < i < 50) 
Fig. 1. Synchronization errors of network (fT6l-(fT7l. 
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Since 

||fii(Xl,X 2 , ■ • ■ ,XiV,S)[|l 

= (/ 1 (e J )-2/ 1 (e J+1 ) + / 1 (e i+2 )) 2 

+ (/3(Xi, Xi+1, X J+2 , S)) 2 

< (o||ei||i + 2a||e i+ i||i + a||e J+2 ||i) 2 

+ (M\\e i \\ 1 + 2M||e i+ i|| 1 + M\\^ +2 \\x) 2 

< 6(a 2 + Af 2 )(|| ei || 2 + ||e s:+1 || 2 + ||e i+2 || 2 ) 

< 18( a 2 + Af 2 )(||e t ||i + ||e 1 . +1 || 2 + ||e 4+2 || 2 ) 

< 18(a 2 + M 2 )(|| ei || 2 + ||e i+1 || 2 + ||e 4+2 || 2 ) 2 , 

where 1 < i < 50, one gets 

||hi (xi,x 2 , • ■ ■ ,X jy,s)|| 2 

< 3V2(a 2 + M 2 )(||e. i || 2 + ||e t+1 || 2 + ||e i+2 || 2 ). 

Thus H2 and H3 hold. According to Theorem 2, the 
synchronous solution S(t) of dynamical network ( I16> - (I17> is 
globally asymptotically stable. 

Assume that k { = 1, d 4 (0) = l,Xj(0) = (4 + 0.5i, 5 + 
0.5i, 6 + 0.5i) for 1 < i < 50 and s(0) = (4, 5, 6). The 
synchronous error e.; is shown in Fig. 1. Obviously, the zero 
error is globally asymptotically stable for dynamical network 

GSJ-O- 

Remark 2: It is well known that the nearest-neighbor cou- 
pled ring lattices are very hard to synchronize. This is because 
the coupling coefficient c satisfies c = 0(N 2 ). However, 
the above example shows that the synchronization of nearest- 
neighbor coupled ring lattice will be relatively easy by adding 
a simple adaptive controller. 
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V. Conclusions 

We have further studied the locally and globally adaptive 
synchronization of an uncertain complex dynamical network. 
Several novel network synchronization criteria have been 
proved by using Lyapunov stability theory. Compared with 
some similar results, our assumptions and adaptive controllers 
are very simple. Furthermore, the effectiveness of these syn- 
chronization criteria have been demonstrated by numerical 
simulations. 
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